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IDCS J1426.5+3508: COSMOLOGICAL IMPLICATIONS OF A MASSIVE, STRONG LENSING CLUSTER ATZ = 1.75 
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ABSTRACT 

The galaxy cluster IDCS J 1426.5+3508 at z — 1.75 is the most massive galaxy cluster yet discovered at z > 1.4 
and the first cluster at this epoch for which the Sunyaev-ZeF Dovich effect has been observed. In this paper we 
report on the discovery with HST imaging of a giant arc associated with this cluster. The curvature of the arc 
suggests that the lensing mass is nearly coincident with the brightest cluster galaxy, and the color is consistent 
with the arc being a star-forming galaxy. We compare the constraint on M200 based upon strong lensing with 
Sunyaev-ZeF Dovich results, finding that the two are consistent if the redshift of the arc is z > 3. Finally, we 
explore the cosmological implications of this system, considering the likelihood of the existence of a strongly 
lensing galaxy cluster at this epoch in a ACDM universe. While the existence of the cluster itself can potentially 
be accomodated if one considers the entire volume covered at this redshift by all current high-redshift cluster 
surveys, the existence of this strongly lensed galaxy greatly exacerbates the long-standing giant arc problem. 
For standard ACDM structure formation and observed background field galaxy counts this lens system should 
not exist. Specifically, there should be no giant arcs in the entire sky as bright in F8 14W as the observed arc for 
clusters at z > 1.75, and only ~ 0.3 as bright in F160W as the observed arc. If we relax the redshift constraint 
to consider all clusters at z > 1 .5, the expected number of giant arcs rises to ~ 15 in F160W, but the number of 
giant arcs of this brightness in F814W remains zero. These arc statistic results are independent of the mass of 
IDCS J1426. 5+3508. We consider possible explanations for this discrepancy. 

Subject headings: galaxies: clusters: IDCS J 1426.5+3508 - gravitational lensing: strong - cosmology: obser- 
vations, cosmological parameters 



1. INTRODUCTION 

Galaxy clusters have historically played a central role in 
cosmology, with the most massive and distant systems pro- 
viding the most profound insights. For example, observa- 
tions of the Coma clu ster provided the first evidence for dark 
matter (IZwickvlll933l) . while the existence of exceptionally 
massive clusters at early times was an important argument 
forOo < 1 (e.g.. ICarlberg et al.lll996t iDonahue et alJll998l 
and references therein). In recent years much attention has 
been given to the question of whether the most distant, highest 
mass clusters are consistent with a standard Gaussian ACDM 
cosmology, or whether one must invoke non-Gaussianity of 
the initial density fluctuations from inflation that seed struc- 
ture for mation. While these analyses have yiel d ed divergent 
results dHovle et all [201 It lEnavist et al.l 120111: ICayon et all 
1201 U [Williams on et al.ll2011l) . it is clear that the most mas- 
sive, distant clusters remain valuable cosmological probes. 

The galaxy cluster that is the focus of this paper was de- 
tected as part of the IRAC Distant Cluster Survey (IDCS), 
an ongoing 8.82 deg 2 s urvey within t he Sp itzer Deep, Wide- 
Field Survey (SDWFS lAshby et all 120091) region that em- 
ploys full photometric redshift probability distributions for a 
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4.5/mi-selected galaxy catalog to identify galaxy clusters at 

< z < 2. This program extends the IRAC Shallow Clus- 
ter Survey (ISCS; Eisenhardt et al. 2008), which has yielded 
the largest sample of spectroscopically confirmed clusters at 

1 < z < 1 .5 by pushing both to lower mass and higher redshift. 
The cluster IDCS J 1426.5+3508 was identified as a strong 
candidate for a high mass, z > 1.5 cluster in this program and 
targeted for detailed follow-up. Spectroscopic observations 
with the HST WFC3 grism a nd LRIS on Keck, described in 
detail in lStanford et al.l d2012l) . confirm that this cluster lies at 
z= 1.75. 

While the existence of a cluster at this redshift is not sur- 
prising, multiple line s of evidence n ow suggest that this is a 
truly massive cluster. Stanfor d et al.l (1201 2b report a Chandra 
X-ray mass estimate of M 200- 5.6 x 10 14 M a , wh ile Sunyaev- 
ZeF Dovich imaging from iBrodwin et all (1201 2l) implies that 
the mass contained within a region overdense by a factor of 
500 relative to critical density is M 50 o= 2.6 ± 0.7 x 10 14 
M G , (M200- 4.2 x 10 14 M Q for a typical halo concentra- 
tion). For comparison, this Sunyaev-ZeF Dovich mass is 
only ~ 40% lower than that of XMMU J22 35. 3-2557 at 
z = 1.39 dMullis et al.ll2005t iRosati et al.ll2009l) . which is the 
only published cluster at z > 1.2 more massive than IDCS 
J1426. 5+3508. Meanwhile, the one spectrosco pically con- 
firme d galaxy cluster at higher redshift (z = 2.07, Gobat et alj 
1201 lb has an estimated total mass of 5.3 - 8 x 10 13 M Q - a 
factor of five to ten lower than IDCS J1426.5+3508. 

In this paper we focus upon the discovery of a giant arc as- 
sociated with this cluster and the implications of its existence 
in the context of ACDM structure formation. The layout of 
this paper is as follows. In Sj2]and[3]we present the discovery 
and attempted spectroscopy of the giant arc. In fj4]we derive 
strong lensing constraints on the cluster mass, and discuss 
the redshift regime over which these constraints are consis- 
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tent with the SZ mass. We then extend our discussion in |5]to 
consider the probability for the discovery of a giant arc associ- 
ated with this cluster. Finally, in f|6]we summarize our results 
and consider potential theoretical modifications that may re- 
solve the arc statistic discrepancy. Throughout this paper we 
use cosmological parameters consistent with the seven year 
WMAP results (Qa = 0.728, H = 70.4 km s -1 , cr 8 = 0.809; 
lKomatsuetal]|2011l) . 

2. DETECTION OF A GIANT ARC 

We obtained HST imaging w ith the Advanced Camera for 
Surveys (ACS: iFord et all 120031) and Wide-Field Camera 3 
(WFC3; iKimble et alj l2008t) as part of Cycle 17 program 
1 1663 (PL Brodwin) between 08 July 2010 and 07 November 
2010. The total exposure times in F814W and F160W were 
4. 5 ks and 2.6 ks , respe ctively. Further details are provided 
in IStanford et al.l (120121) . Subsequent grism observations in 
Cycle 18, coupled with Keck spectr oscopy, confirm that the 
galaxy cluster is real and atz = 1.75 (Stanf ord et al.ll20l2h . 

Within the HST/ACS and WFC3 imaging we iden- 
tify a highly elongated object which we interpret as a 
strong arc lensed by the cluster. We present a composite 
F814W+F160W image of the cluster field in Fig. Q] high- 
lighting this object. The length is 478 but the width is un- 
resolved in the HST photometry, and hence it easily satisfies 
the standard length-to - width criteria for a giant arc (l/w > 10; 
IWu & Ham mer 1993). The curvature of the object is consis- 
tent with lensing by the cluster potential, and the color can be 
used to further constrain the nature of the object. 

2.1. Arc Photometry 

We extract the magnitude of the arc within a polygonal 
aperture constructed to enclose the full extent of the arc. The 
enclosed area within the aperture is 9.0 arcsec 2 , and the ex- 
tent along the major axis is 677 (Fig. ED) . As part of this 
proce dure we first use Source Extractor dBertin & Arnoutsl 
1996) to generate a background map. The flux and back- 
ground are measured both within the source aperture and for 
an ensemble of blank-sky apertures surrounding the arc, from 
which we can calculate the aperture-to-aperture photomet- 
ric variance. We measure F814W= 24.29 + 0.31 mag and 
F160W= 23.75 + 0.21 mag (AB). Given the large uncertain- 
ties in the total magnitude, the integrated color of the arc, 
F814W-F160W= 0.55 ± 0.37 mag, is only modestly con- 
strained. To obtain an improved estimate of the color, we re- 
compute the color within a smaller, 078 x074 rectangular aper- 
ture that includes the region with the highest signal-to-noise 
ratio (Fig. [TJ. Within this aperture we obtain a more precise, 
statistically consistent color, F814W-F160W= 0.25 + 0.13 
mag. As expected for an arc, we see no evidence for color 
changes along its length that might be indicative of a chance 
superposition of sources. 

What does this color imply about the source galaxy? The 
fact that the source is not a drop-out in F814W constrains the 
redshift to be z < 6. Th e observed color is sim ilar to that of 
z ~ 4 B-band dropouts dGonzalez et al.ll201 ll) . For a some- 
what more complete picture of plausible redshifts, we con- 
sider t he predicted colors from a suite of Bmzual & Chariot 
(2003) stellar population models. For a Chabrier (2003) mass 
function and solar metallicity, the observed color is incon- 
sistent with passively evolving stellar populations at all red- 
shifts. To be precise, single burst models in which star for- 
mation has ceased more than -100 Myrs earlier uniformly 
predict F814-F160W colors that are significantly redder than 



observed regardless of the formation redshift. Conversely, 
models with an exponentially declining star-formation rate 
(r-models) can successfully reproduce the observed color at 
higher redshifts. These models place a lower bound on the 
redshift as a function of t and the assumed formation red- 
shift. For t = 1 Gyr, lcr consistency with the observed color 
requires z > 1.75 (2.5) for Zf orm > 2 (3). Smaller values of t 
increase the minimum redshift. 

2.2. Arc Geometry 

The location of the arc is such that no nearby individual 
galaxies appear to be contributing significantly to the lens- 
ing, and the curvature is consistent with the centroid being 
nearly coincident with the brightest cluster galaxy (BCG). If 
we make the assumption that the lensing is indeed centered 
upon the BCG, then the radius of the arc is 6 = 1476 ± 072. 
There is however no guarantee that the BCG lies directly at 
the bottom of the potential in the absence of a detailed lensing 
model, so for the analysis below we assume an uncertainty 
of 30 kpc in the centroid, which corresponds to a 374 uncer- 
tainty in the radius of the arc. For local clusters from the Lo- 
CuSS program roughly 68% of clusters have offsets between 
the BCG and X-ray peak less than this value (Sanderson et al. 
l2009h . 

3. GEMINI AND HST SPECTROSCOPY 

We attempted to obtain a redshift for the giant arc using 
Gemini North. We targeted Lya emission during 6.5 hr of Di- 
rector's Discretionary Time, using GMOS for long slit spec- 
troscopy. The data were acquired between 26 June 201 1 and 
06 July 2011 (UT) using the B600 grism, which is blazed at 
4610 A. We observed with a central wavelength of 5150 A 
and two pixel binning in both the spatial and spectral direc- 
tions using a 170 slit. The resultant spatial and spectral reso- 
lutions are 07146 pix~' and 0.9 A pix~', respectively, and the 
slit was positioned to lie along the long axis of the arc at a 
position angle of 270.85 degrees. We obtained 13 x 1780 s 
exposures, dithered in both the spatial and the spectral dimen- 
sions, for a total on-source exposure time of 6.43 hr. The 
seeing during these observations ranged from 0767 to 1708, 
with a median seeing of 0777 based on an early M-type star 
which was serendipitously observed in the majority of the ob- 
servations. Conditions during the observations were mostly 
photometric, but included some data taken during 70th per- 
centile (patchy cloud) conditions. 

The redshift constraints arising from the cluster redshift 
(i.e., z > 1.75) and the fact that the arc is detected in the 
F814W imaging (i.e., z < 6) provide a first bound on the 
redshift. Within this range we focused upon z < 4 and de- 
signed the observations to be sensitive to strong Lya emis- 
sion (comparable to a Lyman Alpha Emitter), if present for 
2 < z < 4. We reduced the data using standard long slit 
procedures within IRAF. Unfortunately, we detect neither any 
continuum nor any emission lines at 3590 A < A < 6660 A 
(2.0 < z < 4.5 for Lya). 

We also attempted to obtain a redshift for the arc us- 
ing ffi>7YWFC3 grism data from program 12203 (PI Stan- 
fo rd). The data and red uction procedure are both described 
in IStanford et al.l d2012l) . Neither continuum nor emission 
lines for the arc were detected with either the G102 or G141 
grisms, which together cover the wavelength range 0.8/im 
< z < 1 .65/mi. Using simulations run with the aXeSIM soft- 
ware, we calculate that the corresponding emission line de- 
tection threshold corresponds to / < 4 x 10~ 16 erg s _1 cm" 2 
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FIG. 1. — Left- Combined F814W+F160W image of the cluster center and giant arc. The field of view is 30"; North is up and East is to the left. Right- 
Zoomed-in version of the same image centered on the arc. The polygon is the aperture used to extract the arc photometry, while the smaller rectangle is the region 
within which the color was determined. The field of view is 12'.'5 across. The images have been smoothed with 5 pixel and 3 pixel Gaussian kernels, respectively, 
to enhance the contrast. 



at 0.95^m and / < 3 x 1(T 16 erg s -1 cirT 2 at 1.35//m 
(5 <t). 8 Given the non-detection of Lya with either GMOS 
or HST/WFC3 spectroscopy, HST narrow- or medium-band 
imaging may be the most promising avenue for refining the 
redshift estimate for this arc. 

4. CLUSTER MASS FROM STRONG LENSING 

4.1. Enclosed Mass within the Arc 

Under the assumption of circular symmetry for the cluster 
lens, we calculate the total mass enclosed by the giant arc as 
a function of the source redshift. In this case the arc radius f? a 
(~ 125 kpc at zl = 1-75) identifies the radius of the tangential 
critical curve, which can be easily related to the enclosed mass 
through the relation, 

M. d = nZ c ft 2 , (1) 
where E c is the lensing surface critical density, which reads 

c 2 D s 

E c = — . (2) 

47rGD L D LS 

In this equation L\, £>s, and Z\s are the angular diameter 
distances to the lens, to the source, and from the lens to the 
source, respectively. 9 

We emphasize that this enclosed mass is independent of the 
specific density profile assumed for the lens. One important 
caveat in this estimate, however, is that the assumption of cir- 
cular symmetry is known to yield an overestimate of the en- 
closed mass for more realistic systems with intrinsic elliptic- 
ity. To approximately account for this effect, we assume that 
the circular model results in a factor o f ~ 1,6 overestimate of 
the mass, consistent with Bar telmannl (11995), and quote val- 
ues below that include this correction. 

In Figure |2] we show the resulting enclosed mass as a func- 
tion of the source redshift. The closer the source is to the 
deflector, the larger the enclosed mass needs to be due to the 
geometric suppression of the lensing efficiency. The uncer- 
tainty in the enclosed mass shown in the Figure corresponds 

8 |http : / / axe ■ st sci . edu/ axesim/ 

9 We refer the reader to Meylan et al. (2006) for a detailed review of grav- 
itational lensing. 



to the uncertainty in the arc radius, for which we adopt the 
nominal value of 30 kpc (see £|2}- The value for the enclosed 
mass reaches a lower limit of M fl = 6.9 + 0.3 x 10 13 M for 
z s = 6. The enclosed mass in this central 125 kpc region, 
which contains minimal assumptions, already is comparable 
to the total mass inferred for the only spectroscop ic ally con- 
firmed cluster at higher redshift (Gobat et al. 201lT). 

4.2. M 200 

The next step is to estimate the total mass within rajo f° r 
the cluster. This problem is underconstrained, necessitating 
several simplifying assumptions. We initially assume that the 
density profile of t he cluster is well represen ted by a spher- 
ical NFW model dNavarro et al.l fl996l 1 19971) . For a given 
virial mass, we compute the concentratio n of the dark mat - 
ter halo according to the prescription of iGao et all ([2008 ), 
which is in tu rn a modified ver sion of the original NFW pre- 
scription. The lGao et al.l (2008) formula has been extensively 
tested against numerical simulations, including the high red- 
shift regime relevant to the current analysis, and is expected 
to pro vide improvement over the prescriptions of E ke et all 
(l200lh and lBullock et alJ (l200l . 

To account for asymmetries in the cluster mass distribu- 
tion, we next assign a non-vanishing ellipticity to the lens- 
ing potential, accordin g to the procedure summarized in 
Meneghett i et al.l (120031) . Finally, we assume that the arc is 
produced by a source lying near one of the caustic cusps sit- 
uated along the major axis of the lens, so that the arc radius 
corresponds to the maximum elongation of the critical curve. 
Thus, we vary M200 until we find a match between this maxi- 
mum elongation and 6* a . 

In Figure 12 we show the resulting M200 as a function of the 
source redshift. We assume an ellipticity, e m 0.32, consis- 
tent with t he mean of the ellipticity distribution presented in 
Figure 7 of iFedeli & Berciano Albal (120091) . and use the stan- 
dard deviation of this distribution cr e 0.074 to define the 
uncertainty shown by the shaded region. A caveat to this as- 
sumption jsjhjinhisejlipticity distribution is derived at low 
redshift. iLee et al] ([2005) however demonstrated that evolu- 
tion of the ellipticity distribution is expected to be negligible 
forz < 1.5, and sufficiently small at z < 2 as to not impact our 
calculations. From this analysis the derived value for M200 ap- 
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proaches a lower limit of M 2 oo= 2.8^° x 10 14 M as z s — * 6, 
where the quoted uncertainty reflects the uncertainty in the 
ellipticity. For reference, we also show the results obtained 
with the same fiducial ellipticity e m , but assuming the original 
NFW prescription for the concentration (dashed curve). Since 
at high redshift the NFW concentration is always higher than 
other prescriptions, the required M2 00 of the cluster is 3 5% 
smaller. Conver sely, use of either the Bulloc ket al.l (1200 ll) or 
lEke et all ([2001) prescriptions would lead to a larger value of 
M 2 oo. 

It should be noted that strong cluster lenses are usually a bi- 
ased subsample of the whole cluster population, in the sense 
that they tend to be intrinsically more concentrated, and to 
be prolated with the major axis aligned along the line of sight 
dHennawi et alj2007tlMeneghetti et alj2010t) . This bias is ex- 
pected to be even more severe in redshift and/or mass ranges 
where strong lensing is particularly rare, such as the case un- 
der consideration. It has also recently been observed that there 
may be a stronger than expected correlation between con- 
centration and cluster mass, with lower mas s systems having 
higher than expected concentrations (e.g ., ISchmidt & Alien! 
120071; lEttori e"taill20 lot lOguri et alj|2012b . Therefore a con- 
centration higher than that provided by the iGao et all (120081) 
prescription, which is a mean over the entire cluster popu- 
lation, might actually be more realistic in this circumstance. 
As an example, f or a sample of lensing clusters at z ~ 0.5 
lOguri e t al. (2009) found conce ntrations a factor of ~ 2 higher 
than would be predicted by the | Gao et al.l (|2008) distribution. 
Simulations by Meneghet ti et al] (l 2010l) also indicate that pro- 
jection effects may yield observed concentrations that are en- 
hanced by up to a factor of two. In the current case, if we 
assume that the concentration is a factor of two above the 
IGao et ail (120081) prescription, then the derived M200 would 
decrease by roughly a factor of 2.2. 

In Figure |2] and Table Q] we also present the masses and as- 
sociated uncertainties derived from the Sunyaev-Zel'Dovich 
analysis dBrodwin et al .1120121) . The two approaches appear to 
yield consistent estimates for M200 if the source redshift lies 
at z > 3.5, or more conservatively z > 3 if one includes the 
potential for reducing the lensing mass by up to a factor of 
~2 if the halo concentration is larger than for a typical cluster. 
Coupled with the color of the arc, these factors together argue 
that the most plausible redshift is 3 < z < 6. 



5. COSMOLOGICAL IMPLICATIONS 

The redshift of this cluster makes it a unique and interest- 
ing test for cosmological structure formation. To be specific, 
the most distant clusters known to host giant arcs prior to 
this s tudy lie at z ~ 1 (e.g. JGladders et aTll2003b lHuang et al.l 
120091) . IDCS J1426.5+3508 significantly extends the redshift 
baseline over which arcs are known to exist. In this section we 
consider the probability for this massive, strong lensing clus- 
ter to exist and be detected in our survey. Specifically, given a 
standard ACDM cosmology with the seven year WMAP cos- 
mological parameters, what is the probability of detecting a 
giant arc of this brightness behind a cluster at z > 1.75? 

In order to estimate how rare the observed gravitational arc 
is, we evaluate the number of arcs in the whole sky that are 
expected to be produced by galaxy clusters at redshift larger 
than zl = 1-75. We first estimate, for a fixed source redshift 
Zs, the contribution to the optical depth by structures in the 
desired redshift range, 




FIG. 2.— Mass of IDCS J1426.5+3508 as a function of redshift of the 
lensed source. The lower curve corresponds to the mass enclosed within the 
arc, with the shaded region denoting the uncertainty associated with the offset 
of the BCG relative to the clu ster potential. The upper curve is the inferred 
M200 assuming the Gao et al j 1 2008) prescription for the concentration and 
an ellipticity e m ^ 0.32 for the cluster dark matter halo. In this case the 
uncertainty denoted by the shaded region is dominated by the intrinsic scatter 
in the distribution of halo ellipticities. We also overplot as a dashed line the 
inferred mass if one instead uses the original NFW prescription for the halo 
concentration (which can be considered a lower bound). The horizontal line 
and associated uncertainties correspond to the M200 derived from Sunyaev- 
Zel'dovich observations. In this case the uncertainties do not include the 
potential systematic bias associated with extrapolating SZ scaling relations 
to higher redshift. 



Ta(Zs) 



dMn(M,z) 



dV(z) 



dz 



0-a(M, Z ). 



(3) 

In the previous equation n(M, z) is the mass function of cos- 
mic structures, dV(z)/dz represents the comoving volume of 
space per unit redshift, and cr q (M,z) stands for the cross sec- 
tion of individual clusters for images having the morphologi- 
cal property q. In what follows we assume that q is a length- 
to-width ra tio > 10, as customa ry in arc statistics studies, and 
employ the lTinker et al.l (120081) mass function. 

The total number of arcs with the property q that are ob- 
served in the sky with a magnitude brighter than m then sim- 
ply reads 



Nq(m) — An ns(m) 



p(zs,m)T q (zs)dzs, 



(4) 



where p(zs,m) is the source redshift distribution, while ris(m) 
represents the observed number density of sources with mag- 
nitude lower than m, i.e., the cumulative number counts. We 
adopt the redshift distribution a nd number counts for sources 
in the Hubble UDF provided bv lCoe et all j2006h for F775W 
and F160W, and correct the number counts for the lens- 
ing magnification bias using the same procedure deta iled in 
iFedeli et all d2008l) and lFedeli & Berciano Albal (120091) . 

For the practical comp utation of th e optic al depth, we use 
the strategy detailed in Fed eli et ail (|2008) and references 
therein. Briefly, merger trees are constructed based on the 
extended iPress & Schechterld 19741) theory, which represent a 



A Strong Lensing Cluster at z — 1 .75 



5 



100 



I 0.1 - 



0.01 



0.001 




23.5 24 24.5 25 25.5 26 



FIG. 3. — Predicted number of giant arcs over the entire sky as a function 
of magnitude in F775W and F160W for clusters at z > 1.75. The points 
correspond to the results from our simulations, while the curves are spline 
interpolations between the data points. The arrows represent the all-sky lower 
limits derived from the observed arc in IDCS J 1426.5+3508, with the width 
at the bottom of the arrows corresponding to the photometric uncertainty. We 
note that finding one arc per 8.82 deg 2 would correspond to ~ 4700 arcs 
all-sky. 



model of the cluster population. A lensing potential elliptic - 
ity is assi gned to each cluster, ext r acted from the distribution 
shown in Fedeli & Berciano Alba (2009), and cluster dynam- 
ical activity occurring at the knots of the merger trees is suit- 
ably modeled. Individual clusters are modeled as NFW pro- 
files with the assigned potential elliptici ty, and the conc entra- 
tion is linked to the mass through the iGao et all (120081) pre- 
scription for the concentration, consistent with §4.2\ We in- 
clude a lognormal scatte r in the concentration with cr c = 0.2, 
as in iFedeli et"ai1 (120071) . Finally, the cluster cross sections 
for giant arcs are computed using the fast and semi-analytic 
prescription of Fede li et all (12006b . and the optical depth inte- 
grals are approximated by using a Monte-Carlo scheme. 

Figure [3] shows the resulting number of arcs expected 
across the full sky as a function of magnitude in F775W 
and F160W. For our observed arc we assume a color correc- 
tion F775W-F814W- 0.0 (AB), consistent with the expected 
color of a star-forming galaxy at this epoch. For the observed 
magnitudes we expect to find no arcs over the entire sky as 
bright in F814W as the one we observe and only ~ 0.3 as 
bright in F160W. Indeed, none are expected within 0.5 mag- 
nitudes of the brightness of the arc in F814W. Given that the 
area of our survey is only 8.82 deg 2 (2 x 10~ 4 of the full sky), 
the detection of an arc is highly implausible. For reference, 
we test the sensitivity of this result to the redshift of the lens- 
ing cluster. Even if one were to consider the entire cluster 
population at z > 1-5 rather than setting the observed cluster 
redshift as the minimum permitted redshift, then the expected 
number of arcs detected remains zero in F814W and ~ 15 
all-sky in F160W. The latter still leaves only a probability of 
2 x 10~ 3 that we would have detected such an arc in this sur- 
vey. Note that the specific mass of IDCS J1426.5+3508 does 
not enter into our calculations - we posed the question of how 
many arcs should be produced by all clusters. 

At lower redshifts the excess of giant arcs behind clus- 



itraints on the source redshift distrib ution (e.g., 
et al.M2003b IW ambsgans s et al 112004 iDalal et al. 
alJ 120061; IFedeli et all 120081; IWambsganss et al.1 
Disio & Nataraian 2011), the issue is not wholly 



ters relative to pred ictions has been re alized for well over 
a decade (e.g., iBartelmann et alJ 1 19981) and is known as 
the arc statistics problem. While significant effort has 
been devoted to reconcile the observations with improved 
models that incorporate more detailed physics and im- 
proved constraints on t he source redshift distribution (e.g., 
Bartelmann i 
2004 iLi et i 

2008t lr7 Aloisio & Nataraian 2011), the issue is not wholly 
resolved (e.g.. Ilvleneghetti et al.l 120 1 lb . The current cluster 
greatly exacerbates the situation - this arc simply should not 
exist. 

It is interesting to break down the above calculation of the 
number of arcs (Equation [4} to ascertain the dominant factor 
that drives the expected number of arcs to zero. The two fun- 
damental physical quantities that determine the total number 
of arcs are the total lensing optical depth of all clusters and 
the number density of background sources to be lensed. 

The lensing optical depth depends upon the products of the 
cluster mass function and the lensing efficiency of individual 
clusters, which is a strong function of cluster mass. An ob- 
vious explanation would be if IDCS J 1426.5+3508 exceeds 
the mass of clusters expected to exist in a ACDM universe 
at this redshift. If lensing requires reaching a critical mass 
threshold that is exceeded in reality but not in our model, then 
the observed number of arcs would clearly exceed expecta- 
tions. In that case, a plausible solution would be to invoke 
non-Gaussianity to enhance the nu mber of extremely m assive 
clusters. It is however argued in Brodwin e t al.l (1201 2l) that 
IDCS J 1426.5+3508 is consistent with ACDM, indicating that 
invoking non-Gaussianity is not necessarily the appropriate 
solution. 

A more subtle solution would be if the lensing cross- 
sections of individual clusters systematically exceed the cal- 
culated theoretical cross-sections. There is evidence that 
this lensing efficiency is indeed systematically underesti- 
mated. Several papers have argued that including the im- 
pact of baryonic contraction can raise the lensing cross- 
section by between 25% a nd a factor of two dLi et al.ll2006t 
IWambsganss et al.ll2008h . iMeneghetti et all (1201 lh also find 
that simulated clusters produce ~ 50% fewer arcs than X-ray 
clusters at z ~ 0.5 - 0.7. While these corrections work in 
the proper direction, it appears that even if one imposes rea- 
sonable physical tweaks to account for the impact of baryonic 
infall on the density profile, the change remains insufficient 
to account for th e arc in IPCS J14 26. 5+3508. Moreover, as 
demonstrated by iMead et alJ J2010). inclusion of AGN feed- 
back acts to counteract the impact of baryonic infall, resulting 
in a smaller enhancement to the cross-section. Even ignoring 
the mitigating effect of feedback, the discrepancy is simply 
too large. Specifically, doubling the predicted number of arcs 
would imply ~ 1 arc all-sky with the observed F160W mag- 
nitude, and would still imply none with the observed F775W 
magnitude. 

Another means of boosting the effective cross-section for 
an individual system is via the presence of addi t ional struc- 
tures along the line of sight. Puch weiii & Hilbertl (2009) used 
the Millenium simulations to quantify the impact of such sec- 
ondary structures. These authors found that the typical impact 
is to enhance the cross section by 10 - 25%, with enhance- 
ments of 50% not uncommon for individual systems. Again, 
this factor alone is insufficient to explain the discrepancy. 

The final factor that can drive the prediction of zero arcs is 
underestimation of the surface density of galaxies sufficiently 
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Table 1. Derived Masses for IDCS J1426. 5+3508 



Method 


Radius 


Mass (10 14 Ms) 


Lensing 


125 kpc 


> 0.69 ± 0.03 


Lensing 


''200 




SZ 


'500 


2.6 ± 0.7 


sz 


loo 


4.1 ± 1.1 



NOTE. — The lensing limiting masses correspond to the maximum pos- 
sible source redshi ft, z s = 6 - Both M200 masses presume a concentration 
consistent with the Gao et al. (2008) relation. If the projected concentration 
is a factor of two higher then the lensing and SZ masses drop by factors of 2. 1 
and 1.2, respectively. As noted in the text, r^QQ and r2oo are defined relative 
to critical density. 



bright to yield a source similar to the observed arc after mag- 
nification by the cluster potential. One concern here is that 
we have used results from the UDF to inform our redshift 
distribution for background sources, yet this region is suffi- 
ciently small that cosmic variance is a concern. As an alter- 
nate test, we also try an analytic prescription for the back- 
ground nurribe£cojmtsbased upon the z ~ 3 luminosity func- 
tion of iReddv et alj (120081) to model the background distri- 
bution. We find however that this approach does not quali- 
tatively alter the results of our analysis. It therefore seems 
unlikely that error in the source distribution is the origin of 
the discrepancy. 

We therefore identify no obvious physical solution to ex- 
plain the existence of this arc, though exceptionally high con- 
centrations seems like the most promising avenue to explore. 
Finally, a last possible solution would be if the observed 
source is not a background arc. Given the combination of 
color, curvature, l/w, and lack of color variation, this possi- 
bility also seems unlikely. 

6. CONCLUSIONS 

We have presented evidence for the existence of a giant 
arc behind the massive galaxy cluster IDCS J 1426.5+3508 at 
Z = 1.75. This unique system constitutes the highest redshift 
cluster known to host a giant arc. From the strong lensing we 
derive an enclosed lensing mass within the central 125 kpc of 
M. d > 6.9±O.3xlO 13 M , which is comparable to (or exceeds) 
the total masses of most other known clusters at z > 1.5, 
confirming that IDCS J 1426.5+3508 is an exceptionally mas- 
sive cluster at this epoch. Having derived the enclosed mass, 
we next provide a comparis on with the Sunyaev-Z erDovich- 
derived value of M200 from iBrodwin et all (1201 2b . We find 
that the two are consistent if the arc lies at z > 3, and derive a 
lower bound on the mass of M 2 oo > 2.8 x 10 14 M . 

Finally, we investigate the cosmological implications of this 
system. The greatest challenge posed by this cluster is ex- 
plaining the existence of the giant arc at all. In $5]we demon- 
strate that under realistic assumptions for the lensing cross- 
section, cluster mass function, and background galaxy dis- 
tribution, the total number of giant arcs behind clusters at 



z > 1-75 that are at least as bright as the observed arc is zero 
in F814W and < 1 in F160W. Very simply, the arc we have 
discovered behind IDCS J 1426.5+3508 is not predicted to ex- 
ist. If one considers an ensemble of lensing clusters extending 
to somewhat lower redshift, z > 1.5, the tension with theoret- 
ical expectations decreases slightly forF160W, but in F814W 
the number of predicted arcs remains zero. We briefly discuss 
possible explanations for this discrepancy, but find no obvious 
solution. A tendency of higher concentrations in observed 
clusters than simulated systems has the greatest potential to 
decrease the disparity, but is unlikely to be sufficient to recon- 
cile theoretical models for arc statistics with the existence of 
a lensing cluster at z = 1.75. 

Looking towards the future, a statistical sample of the most 
massive clusters in the Universe at z = 1-2 will provide 
the means with which to ascertain the true frequency of arcs 
behind high-redshift clusters. Much as the frequency of strong 
lensing clusters at z — 0.5 - 1 was a surprise a decade ago, 
it appears that this higher redshift regime is poised to yield 
further unexpected discoveries. 
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